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Abstract 

Using a 492 fb -1 data sample collected near the T(4S*) resonance with the Belle detector at 
the KEKB asymmetric-energy e + e~ collider, we observe the decay B° — > ppK*° with a branching 
fraction of (1.181q 25 (stat.) ± 0.11(sysi.)) x 10~ 6 . The statistical significance is 7.2<r for the signal 
in the low mass region. We study the decay dynamics of B° — ► ppK*° and compare with f? + — > 
ppK* + . The if* meson is found to be almost 100% polarized (with a fraction of (101 ± 13 ± 3)% 
in the helicity zero state), while the K* + meson has a (32 ± 17 ± 9)% fraction in the helicity zero 
state. The direct CP asymmetries for B° — > ppK*° and i? + — ► ppK* + are measured to be — 0.08± 
0.20 ±0.02 and —0.01 ± 0.19 ± 0.02, respectively. We also study the characteristics of the low mass 
pp enhancements near threshold and the associated angular distributions. In addition, we report 
improved measurements of the branching fractions B{B + — ► ppK* + ) = (3.381q'qq ± 0.39) x 10~ 6 
and B(B° — > ppK°) = (2.51^0 29 =t 0.21) x 10~ 6 , which supersede our previous measurements. 
PACS: 13.25.Hw 

PACS numbers: 
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After the first observation of the charmless baryonic B meson decay, B+ -> ppK+ [l|,l2j, 
many three-body charmless baryonic decays were found 0, 0, 0, 0]. One important and 
intriguing feature of these decays is that the baryon-antibaryon mass distributions all peak 
near threshold. However, the BaBar collaboration recently reported evidence of the decay 
B° — > ppK* but could not establish either the presence or absence of such a threshold 
enhancement [3]. On the theoretical side, it is generally believed that the B — > ppK* decays 
proceed predominantly through a b — > s penguin loop diagram, which could be sensitive to 
new physics from heavy virtual particles in the loop. Large direct CP violation, ~ 20%, 
is predicted using an effective-amplitude approach in the standard model 0. From a pole 
model 0], it is expected that B(B + — > ppK* + ) < B(B + — ► ppK + ) due to the absence of 
some QCD penguin and electroweak penguin contributions in the ppK* + mode, and that 
B(B° — ► ppK* ) < B(B + — ► ppK* + ) due to the absence of a specific pole contribution and 
the external emission diagram in the ppK* mode. 

In this paper, we study the three-body charmless baryonic decays B° — ► ppK*°(K*° — > 
-fT + 7r~) and 5 + — > ppK* + (K* + — > f^7r + ). The polarization of the if* meson is deter- 
mined, which provides information about the relative importance of penguin and external 
W-emission contributions [10]. The differential branching fractions as a function of the 
baryon-antibaryon mass and the polar angle distributions of the proton in the baryon- 
antibaryon system are also presented. The direct CP violation parameters of these two 
decays are also measured. We use a 492 fb _1 data sample, consisting of 535 xl0 6 BB pairs, 
collected with the Belle detector at the KEKB asymmetric-energy e + e~ (3.5 on 8 GeV) 
collider The Belle detector is a large-solid-angle magnetic spectrometer that consists of 
a silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel 
threshold Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation 
counters (TOF), and an electromagnetic calorimeter composed of CsI(Tl) crystals located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons. 
The detector is described in detail elsewhere (l2| . 

The event selection criteria for the primary charged tracks can be found in Ref. (l3| . 
Kg candidates are reconstructed as 7r + 7r~ pairs with an invariant mass in the range 490 
MeV/c 2 < M w + W - < 510 MeV/c 2 . The candidate must have a displaced vertex and flight 
direction consistent with a Kg originating from the interaction point. We use the selected 
kaons and pions to form K* + (— > K s n + ) and K*° (— > K + ir~) candidates. Events with 
a K* candidate mass between 0.6 GeV/c 2 and 1.2 GeV/c 2 are used for further analysis. 
Candidate B mesons are reconstructed in the B° — > ppK* and B + — > ppK* + modes. We 
use two kinematic variables in the center-of-mass (CM) frame to identify the reconstructed 
B meson candidates: the beam energy constrained mass Mt, c = E^ eam — p\, and the 
energy difference AE = Eb — E^ cam , where -Ebeam is the beam energy, and ps and Eb 
are the momentum and energy, respectively, of the reconstructed B meson. The candidate 
region is defined as 5.2 GeV/c 2 < M bc < 5.3 GeV/c 2 and -0.1 GeV < AE < 0.3 GeV. 
The lower bound in AE is chosen to exclude possible background from baryonic B decays 
with higher multiplicities. From a GEANT [14| based Monte Carlo (MC) simulation, the 
signal peaks in a signal box defined by the requirements 5.27 GeV/c 2 < M^ c < 5.29 GeV/c 2 
and \AE\ < 0.05 GeV. To ensure the decay process is genuinely charmless, we apply charm 
vetoes. The regions 2.850 GeV/c 2 < M pP < 3.128 GeV/c 2 and 3.315 GeV/c 2 < M pP < 3.735 
GeV/c 2 are excluded to remove background from modes with r] c ,J/ip and ip',Xco,Xci,h c 
mesons, respectively. The region 2.262 GeV/c 2 < M pK o , M pK - n + < 2.310 GeV/c 2 is also 
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excluded to remove a possible A+ background. From a study of a charmless B decay MC 
sample, there are non-negligible backgrounds in the candidate region due to B + — > ppK + 
and B° — > ppKg. We remove the B candidates when their Mb c and AE values reconstructed 
for the ppK hypothesis are in the signal box. 

After the above selection cuts, the background in the fit region arises dominantly from 
continuum e + e~ — > qq (q = u, d, s, c) processes. We suppress the jet-like continuum 
background relative to the more spherical BB signal using a Fisher discriminant [l5| that 



combines seven event shape variables, as described in Ref. 16]. Probability density func- 
tions (PDFs) for the Fisher discriminant and the cosine of the angle between the B flight 
direction and the beam direction in the T(4S) rest frame are combined to form the signal 
(background) likelihood L s (£&). The signal PDFs are determined using signal MC simula- 
tion; the background PDFs are obtained from the sideband data: 5.23 GeV/c 2 < M^ c < 5.26 
GeV/c 2 and \AE\ < 0.06 GeV for the ppK* mode; 5.25 GeV/c 2 < M bc < 5.26 GeV/c 2 and 
\AE\ < 0.2 GeV for the ppK* + mode. The different selections for sideband regions of the 
two K* modes ensure similar statistics to determine the background PDFs. We require the 
likelihood ratio TZ = C S /(C S + £5) to be greater than 0.7 for both decay modes. These 
selection criteria are determined by optimization of n s /y/n s + rib, where n s and ri^ denote 
the expected numbers of signal and background events in the signal box, respectively. We 
use the branching fractions from our previous measurements |4J in the calculation of n s and 
use the number of sideband events to estimate n^. If there are multiple B candidates in a 
single event, we select the one with the best \ 2 value from the vertex fit. The fractions of 
events that have multiple B candidates are 21% and 32% for the ppK* and ppK* + modes, 
respectively. 



We perform an unbinned extended likelihood fit that maximizes the likelihood function 



-{n K *+n K K+n,qq) - x 



Y[ (n K *P K * + n Kn P K7I + riqqPqg 



M 

1=1 

to estimate the signal yield of ppK* in the region —0.1 GeV < AE < 0.3 GeV, 5.2 GeV/c 2 < 
M bc < 5.3 GeV/c 2 and 0.6 GeV/c 2 < M Kn < 1.2 GeV/c 2 ; here N is the number of events 
in the fit, and riK*,nKn and n qq - are fit parameters representing the yields of B — > ppK*, 
B — > ppKir and continuum background, respectively. Each PDF is the product of shapes in 
Mb c , AE and M^, which are assumed to be uncorrelated, e.g. for the ith event, P p pK* = 
P Mhc (M bCi ) x P AE {AEj) x P Kw (M KTt ) . 

For the PDFs of ppK* , ppK* + and ppKir decay modes, we use a Gaussian function to 
represent PM hc and a double Gaussian for Pae with parameters determined by MC signal 
events. Moreover, we use a p-wave Breit-Wigner function [TtJ to parameterize the Pm Kit dis- 
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tribution for ppK and ppK* + and use a function obtained by the LASS collaboration [18 
for ppKir. The parameters of these PDFs have been modified to account for the differences 
between data and MC using control samples of J/i/jK*° and J/ipK* + with J/ip pp. The 
modifications related to the mass peaks are all less than 1 MeV/c 2 . The AE distribution 
has a ~ —3 MeV shift while the modification for its width is ~ 1 MeV. For the continuum 
background PDFs, we use a parameterization that was first employed by the ARGUS collab- 



oration 19j, /(Mbc) oc MbcVl — x 2 e^^ 1 ^ x \ to model the Pjw b0 with x given by M bc /-Ebeam 
and where £ is a fit parameter. The Pae distribution is modeled by a normalized second- 
order polynomial whose coefficients are fit parameters. The PDF Pm K7T is modeled by a 
p-wave function and a threshold function, Pm k „ = r x P p _ wavc + (1 — r) x P t hreshoid and 

Phrcshold OC (M Kn -M K - M n ) s X e [cix(M Kn -M K -M^)+c 2 x(M Kn -M K -M^] where &nd 
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FIG. 1: Distributions of AE (with 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 and 0.812 GeV/c 2 < 
Mk-tt < 0.972 GeV/c 2 ), M bc (with |A£?| < 0.05 GeV and 0.812 GeV/c 2 < M Kn < 0.972 GeV/c 2 ) 
and M KlT (with |A£| < 0.05 GeV and 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 ), respectively, with 
proton-antiproton pair mass less than 2.85 GeV/c 2 for (a) ppK* and (b) ppK* + modes. The 
solid curves, solid peaks, dotted curves and dashed curves represent the combined fit result, fitted 
B — » ppK* signal, B — ► ppKir signal and fitted background, respectively. The areas of dotted 
curves are about 15% of those of the solid peaks. 



C2 are fit parameters. Figure [T] shows the fits used to obtain the B — > ppK* yields in the 
proton-antiproton mass region below 2.85 GeV/c 2 , which we refer to as the threshold-mass- 
enhanced region. The signal yields are 70.li 13 J and 54.2ljgj with statistical significances of 
7.2 and 8.8 standard deviations for the ppK* and ppK* + modes, respectively. The signifi- 
cance is denned as a/— 21n(Lo / L mm ) , where Lq and L max are the likelihood values returned 
by the fit with the signal yield fixed to zero and at its best fit value. 

We determine the angular distribution of the K* meson in the region M p p < 2.85 GeV/c 2 
using likelihood fits to obtain signal yields in bins of cos 9k, where 9k is the polar angle of 
the K meson in the K* helicity frame. The theoretical PDF for the K* meson has the form 
3/2 cos 2 9k for a pure helicity zero state and 3/4 sin 2 9k for a pure helicity one (±1) state. 
We use MC simulation to obtain the efficiency function and convolve it with the theoretical 
forms in order to obtain the final PDFs for different helicity states. The signal yields in bins 
of cos 9k are then fitted with the above two different PDFs where the fraction of the helicity 
zero state is floated in the fit and the total yield is fixed to the experimental result. The B 
yield distributions in bins of cos 9k with the corresponding fit curves are shown in Fig [21 We 
find that the K*° meson has a fraction of (101 ± 13 ± 3)% in the helicity zero state and the 
K* + meson has a (32 ± 17 ± 9)% fraction in the helicity zero state. It is interesting to note 
that the helicity zero amplitude is expected to be dominant in the b —>■ s penguin transition 
due to the (V — A) nature of the weak interaction and helicity conservation in the strong 
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FIG. 2: B yield distributions as functions of cos 9k with fit curves overlaid for (a) the ppK*° mode 
and (b) the ppK* + mode. The fraction of the signal in the helicity zero state is the fit parameter 
and is denoted by Hq. The asymmetries in the fit curves are due to detection efficiencies. The 
underlying theoretical distributions are symmetric. 



interaction [lOj. The systematic uncertainty is obtained from the B — > J/ipK* } J/ip — > 
control sample. We compare our measured K* polarization in the helicity zero state with 
the PDG value [20]. The difference is added in quadrature with the PDG error and the fit 
error to extract the final systematic uncertainty. These uncertainties are 0.03 and 0.09 for 
the ppK* and ppK* + modes, respectively. 

Since the detection efficiency depends on M p p, we determine the B — > ppK* yields in bins 
of M p p. We generate a large phase-space MC sample in order to estimate the efficiencies 
properly where the sub-decay branching fractions of K* to corresponding final states are 
included. The Kir angular distribution is fixed by the measured K* polarization for all 
Mpp bins. The partial branching fractions are obtained by correcting the fitted B yields 
for the mass- dependent efficiencies. The differential branching fractions as a function of 
the proton-antiproton mass for both ppK* and ppK* + modes are shown in Fig. [3], and 
the measured branching fractions for different Mpp bins are listed in Table [H Applying 
3.075 GeV/c 2 < M pP < 3.117 GeV/c 2 for J/ip selection, we find good agreement, within 



la - , between our branching fraction measurements and the PDG values [201 ] . In contrast 
to Ref. Q, we find that a threshold enhancement is present for the B° — > ppK* decay. 
With the charmonium regions excluded, we sum these partial branching fractions to obtain: 
B(B° -> ppK* ) = (1.18±g;|§=b0.11) x 10" 6 and B(B + -> ppK*+) = (3.38±g;^±0.39) x 10~ 6 . 
As a by-product of our analysis, we also use the B + — > ppK + and B° — > ppK$ signals to 
estimate the corresponding branching fractions in different M p p bins. The total charmless 
branching fraction B(B + — > ppK+) is (5.36i '22) x 10~ 6 , which agrees well with our latest 
results, (5.54;°; 2 £ ± 0.36) x 10" 6 [l|. The measured value of B(B° -»• ppK ) is (2.51±g;|| ± 
0.21) x 10~ 6 . This result also supersedes our previous measurement |4|. With improved 
experimental accuracy, the following relationships B{B + — > ppK + ) > B(B + — > ppK* + ) and 
B{B + — > ppK* + ) > B(B° — > ppK* ) are established. These inequalities agree with the 
pole model predictions [9j, but the measured B(B° — > ppK* ) is about a factor of 20 larger 
than predicted. This may indicate that the relative weights of different pole contributions 
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FIG. 3: Differential branching fractions for (a) the ppK*° and (b) the ppK* + modes as a function 
of proton-antiproton invariant mass. Note that the two shaded mass bins contain charmonium 
events and are excluded from the charmless signal yields. The data points for the 2.85 GeV/c 2 < 
Mpp < 3.128 GeV/c 2 mass region are off-scale. 



in Ref. [9j are incorrect. 

Systematic uncertainties are determined using high-statistics control data samples. For 
proton identification, we use a A — ► pn^ sample, while for K/tt identification we use a 
D* + — > D°7i + , D° — > K~7i + sample. Note that the average efficiency difference for PID be- 
tween data and MC has been corrected to obtain the final branching fraction measurements. 
The corrections are about 11.5% and 11.7% for the ppK*° and ppK* + modes, respectively. 
The uncertainties associated with the PID corrections are estimated to be 4% for two pro- 
tons and 1% for one kaon/pion. The tracking uncertainty is determined with fully and 
partially reconstructed D* samples. It is about 1% per charged track. The uncertainty in 
Kg reconstruction is determined to be 4% from a sample of D~ — > K s n~ events. The TZ 
continuum suppression uncertainty of 2.3% is estimated from control samples with similar 
final states, B — > J/ipK* with J/ip — > /i + pT . The uncertainties in the best B candidate se- 
lection are estimated to be 2.0% and 3.5% for the ppK*° and ppK* + modes, respectively, by 
taking a difference in the branching fractions with and without the best candidate selection. 
A systematic uncertainty of 5.2% in the fit yield is determined by varying the parameters 
(or changing the functional forms) of the signal and background PDFs. The MC statistical 
uncertainty is less than 3%. The efficiency error caused by the K* polarization modeling is 
estimated to be 2.4% and 4.0% for the ppK*° and ppK* + modes, respectively, by changing 
the polarization value by ±la. The error on the number of BB pairs is 1.3%, where we 
assume that the branching fractions of T(4S) to neutral and charged BB pairs are equal. 
We first sum the correlated errors linearly and then combine them with the uncorrelated 
ones in quadrature. The total systematic uncertainties are 9.7% and 11.6% for the ppK*° 
and ppK* + modes, respectively. 
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TABLE I: Signal yields and branching fractions B 
ppK*° (left) and B + -> (right). 



(1CT 6 ) in different M p p regions for B 



M pp (GeV) 




Yield B (10- 6 ) 



< 2.0 

2.0-2.2 

2.2 - 2.4 

2.4-2.6 

2.6-2.85 

2.85 - 3.128(veto) 

3.128 - 3.315 

3.315 - 3.735(veto) 

> 3.735 



+4.4 
3.7 

+7.1 
-6.3 
+5.4 
4.5 

+3.3 
2.5 

+4.8 
3.9 

+9.8 
-9.0 
-2.1 

-1.5 



7 1+4.8 
'•1-4.1 



9.0 
25.1 
6.4 
4.5 
9.6 
55.7 
1.5 



2.5 



+2.9 
-2.0 



+0.21 
0.18 

+0.36 
-0.32 

+0.31 
0.26 

+0.22 
0.17 

+0.31 
0.25 

+0.65 
-0.59 

-0.15 
-0.11 



0.58l°;|° 

-0.31 
-0.22 



0.43 
1.28 
0.37 
0.30 
0.62 
3.66 
0.11 



0.27 



Charmless 



7^ 4+17-1 
/0 - 4 -14.7 



i 1O+0.29 



3.38; 



-0.73 
-0.60 



We study the proton angular distribution in the proton-antiproton helicity frame with 
M p p < 2.85 GeV/c 2 . The angle 6 P is defined as the angle between the proton direction and 
the positive strangeness K* (i.e. K*~ or K*°) direction in the proton-antiproton pair rest 
frame. The cos 8 P distributions, shown in Fig. HJ do not have a prominent peaking feature 
toward cos8 p ~ 1, which was first observed in the decay B + — > ppK + j2l|. However, current 
statistics are inadequate to draw any definitive conclusions about B — > ppK*. 

To examine the prediction [8[ that direct CP violation in B + — > ppK* + can be as large 
as ~ 20%, we define the charge asymmetry A c h as (iVj, — Ni)/(Nb + A5) for the ppK* and 
ppK* + modes, where A" and b stand for the efficiency corrected B yield and quark flavor, 
respectively. The results are -0.08 ± 0.20 ± 0.02 and -0.01 ± 0.19 ± 0.02 for the ppK* 
and ppK* + modes, respectively. The systematic uncertainty is estimated from the measured 
charge asymmetry for the sideband data. 

In summary, using 535 xl0 6 BB events and applying charmonium vetoes, we observe the 
B° -> ppK* decay with a branching fraction of (l.lStoil{stat.) ± 0.11(syst.)) x 10" 6 . The 
signal yield is 7O.IIJ3J with a significance of 7.2 standard deviations in the M pP < 2.85 
GeV/c 2 mass region. The K*° meson is found to be (101 ± 13 ± 3)% in the helicity zero 
state, compared to (32 ± 17 ± 9)% for the K* + meson. The smaller K* + polarization in the 
ppK* + decay may be attributed to an additional contribution from external W emission. 
We also observe a low mass pp enhancement near threshold for the ppK* mode. The direct 
CP asymmetries for ppK* and ppK* + are measured to be —0.08 ± 0.20 ± 0.02 and —0.01 ± 
0.19 ± 0.02, respectively. With improved experimental accuracy, the relationships B(B + — > 
ppK + ) > B{B+ ppK* + ) and B(B + -> ppK* + ) > B(B° ppK* ) are established. 
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FIG. 4: Distributions of efficiency corrected signal yields vs. cos 9 P in the proton-antiproton system 
with M pP < 2.85 GeV/c 2 for (a) B° -► ppK*° and (b) B+ -► 
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